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We explore the challenges posed by the violation of Bell-like inequalities by rf-dimensional systems exposed 
to imperfect state-preparation and measurement settings. We address, in particular, the limit of high-dimensional 
systems, naturally arising when exploring the quantum-to-classical transition. We show that, although suitable 
Bell inequalities can be violated, in principle, for any dimension of given subsystems, it is in practice increas¬ 
ingly challenging to detect such violations, even if the system is prepared in a maximally entangled state. We 
characterize the effects of random perturbations on the state or on the measurement settings, also quantifying 
the efforts needed to certify the possible violations in case of complete ignorance on the system state at hand. 

PACS numbers: 03.65.Ud, 03.67.Mn 


I. INTRODUCTION 

Year 2014 has marked the 50 th anniversary of John S. Bell’s 
paper JT) on the Einstein-Podolski-Rosen paradox 0, a land¬ 
mark achievement that has shaped the research in modern 
quantum mechanics. Since then, Bell inequalities, which have 
proven the incompatibility of quantum mechanics with any lo¬ 
cal realistic theory, have been established as a viable tool for 
deepening our understanding of the nature of reality and, from 
a more pragmatic perspective, identify unambiguously corre¬ 
lations of a non-classical nature 0. 

Non-locality tests for discrete-variable systems have been 
formulated, over the years, based on both dichotomic and 
multichotomic local measurement settings 0, showing that 
bipartite cZ-level entangled states are in conflict with local re¬ 
alistic assumptions more strongly than entangled qubits. In¬ 
equalities for bipartite cZ-dimensional systems have been put 
forward 0, while various formulation for multipartite d = 2 
systems have been proposed 0, and extended to generic arbi¬ 
trary dimensional systems (7). Continuous-variable Bell-like 
inequalities have been derived for bipartite states using the 
phase-space formalism |8), and extended to multipartite sce¬ 
narios recently 0. 

The investigation on non-locality of discrete larg e-d sys¬ 
tems is particularly relevant in light of the goal of inferring 
non-classical effects in increasingly large systems. The pos¬ 
sibility to unambiguously signal non-classical correlations, 
strong enough to violate a Bell-like inequality, would pro¬ 
vide undoubtable evidence of the relevance of quantum the¬ 
ory in the large-scale domain, and bring new perspectives to 
the investigation of the quantum-to-classical transition DSD- 
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First steps along similar lines of investigation have been made 
by addressing the class of states dubbed as micro-macro 
ones CD, which have then been extended to macro-macro 
scenarios mi. 

When tackling the phenomenology of non-classical corre¬ 
lations in large-scale systems, the detrimental effects of im¬ 
proper measurements should be considered on the same foot¬ 
ing as non-ideal state preparation. It has been shown that 
coarse-grained measurements do indeed affect the chances of 
falsifying local realistic models, even when dealing with oth¬ 
erwise ideal state resources mi. On the other hand, a typ¬ 
ical assessment of the robustness of Bell-like tests is based 
on the amount of white noise that a given resource state is 
able to tolerate before losing any discrepancy with local hid¬ 
den variable (LHV) models 0. This approach was the basis 
for the conclusion drawn by Kaszilowski et al. on bipartite d- 
level systems 0- However, while such analyses provide use¬ 
ful information on the effects that imperfections and environ¬ 
mental mechanisms have on the revelation of non-classicality, 
a systematic study of the scaling behavior of non-idealities 
against the size of a given system has only been partially ad¬ 
dressed EMU- Yet, this is key information for the quest 
of observing non-classicality at the large-scale limit. For in¬ 
stance, as we address in this paper. Bell inequalities formu¬ 
lated for larg e-d systems require a large number of measure¬ 
ment settings: their inadequate arrangement could have rather 
significant implications for the successful falsification of LHV 
models even when states endowed with a large degree of quan¬ 
tum correlations are considered. In this sense, our study goes 
well beyond currently available literature 0. 

Our study here addresses precisely this important issue, 
which we tackle in a careful quantitative manner. By address¬ 
ing experimentally relevant imperfections in the Bell test of 
d-dimensional systems, we point out, through a careful sta¬ 
tistical analysis, the existence of thresholds in the amount of 
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randomness that a given resource state (or, equivalently, mea¬ 
surement setting) can tolerate. As we find out, such threshold, 
which scales unfavorably with the dimensions of the local sys¬ 
tems being used, is rather severe and poses significant ques¬ 
tions on the suitability of Bell-like tests in systems verging to 
the mesoscopic or macroscopic scale, when evaluated against 
the actual value taken by d. In some conditions, even a max¬ 
imally entangled two-qutrit state could fail to violate a Bell 
inequality built from slightly improper measurement settings. 

The remainder of this paper is organized as follows. In 
Sec. [II] we introduce the inequalities addressed in this work, 
setting the notation that will be used throughout our analy¬ 
sis. Sec. [Ill] is devoted to the analysis of the robustness of 
such Bell-like tests for systems of increasing dimensionality 
and exposed to random perturbations of the otherwise ideal 
measurement settings required by each of the inequalities that 
we study. In Sec.[lV]we extend our study to the case of ran¬ 
dom initial states, showing the inherent difficulties related to 
running Bell tests on such resources, and discussing the effec¬ 
tiveness of more sophisticated statistical sampling techniques 
for ascertaining non-locality in randomly picked states of d- 
dimensional systems. Finally, in Sec. [V] we draw our conclu¬ 
sions. 


II. GENERALIZED BELL INEQUALITIES 


and similarly 


d -1 


P(Ai =Bj + k) = E P{Ai = £,Bj = £ + k mod d). 

e=o 


(3) 


In general, one can show that for local theories only three of 
the four equalities in (|2| can be satisfied at the same time, that 
is for local theories I < 3, while non-local theories could reach 
a value of I = 4. 

Beside the generalization of the standard Bell inequality 
presented above, we also consider a second inequality, which 
at variance with the previous one is built with a number of 
terms that grows with the local dimension of the subsystems. 
Explicitly 


W 21-1 / 2k 

id= E (i- A 

k =0 


d — 1 

P(Ai =Bi+k)+P(Bi =A 2 + k+ 1) 
-\-P(A 2 = B 2 +k) +P(B 2 =A\ k) 

P(A\ = B l -k-l)+P{B 1 =A 2 -k ) 

+P(A 2 =B 2 -k~l)+P{B 2 =A 1 ~k-l] 


■ (4) 


In this Section we review the formulation of Bell inequal¬ 
ities for (/-dimensional systems that are used hereafter. The 
original formulation of the Bell inequalities relates the joint 
probabilities for local measurements with two possible out¬ 
comes (d = 2), typically the polarization of photons or a spin- 
1/2 system (dimension d and the spin quantum number / obey 
the relation d = 21 + 1). However, it is possible to extend (in 
a non unique way) such formulation to systems with higher 
local dimension 0: We consider a coupled bipartite system 
and two parties (Alice and Bob) each of them having access 
to a single partition, of dimension d, of the system (either sub¬ 
system PL or ®) and able to measure two independent observ¬ 
ables, Ai, A 2 and B\, B 2 , respectively. The operators As and 
Bs are assumed to have the same spectrum, A,)Af) = £\Af) 
and Bj\B l j) = £\Bf), with £ = 0,... ,d — 1. Clearly, for d- 
dimensional subsystems the measurement can have d differ¬ 
ent outcomes. The generalized Bell inequalities are obtained 
summing the probabilities that Alice and Bob obtain equal (or 
correlated) outcomes of their measurements. In particular, the 
first Bell expression we will consider reads as 

1 =P(A\ = B\) +P(B{ = A 2 + 1) (1) 

+ P{A 2 = B 2 ) +P(B 2 = Ai), 

where P(Aj = Bj) are the probability for Alice and Bob to 
have the same measure outcomes: 

d~ l 

P(Ai = Bj) = E P(Ai = £,Bj = £), (2) 

e=o 


As before, one can show that the maximal possible value 
achievable by this quantity is Id = 4. However, LHV mod¬ 
els allow only for Id < 2. 

Generalized maximally entangled (Bell) quantum states of 
two (/-dimensional systems, 

\® + ) = ^ Cl t\jh®\j)B, ( 5 ) 

va j =o 

violate these inequalities for any value of d 0- It has been 
shown a that, for | < t> + ), the observables that are suited to test 
the generalized Bell inequalities have associated eigenstates 
of the form 

\k) A ,a = I ex P [if j(k+ (*„)] | j) A , (6a) 

v j=o 

\l)B,b £ ex P [ifj(~ l + $b)\ I J)b- (6b) 

v j=o 

Here, a,b = 1,2 with oti = 0, a 2 = 1/2, p! = ~p 2 = 1/4 0. 
Projections of the state of the system onto such states always 
leads to a violation of I and Id, independently of the particular 
value of d. This is shown in Fig. |T| where we report the exact 
numerical values for both I and /,/ as a function of the spin 
quantum number l of the system. Note that, if the probabilities 
are computed according to the rules of quantum mechanics, 
the maximal value of I achievable for d = 2 is / max =2+ \/2 « 
3.4i 1321 and that numerical data show that I monotonically 
decreases with d. 
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FIG. 1. (Color online) Values for the generalized Bell expressions 
I and plotted against the spin quantum number / = (d — l)/2 of 
each subsystem. In order to evaluate I and Id we have used the 
state I®" 1 ") and the measurement operators given in Eqs. |6j. Bot¬ 
tom (green) lines: Maximum values achievable through local hidden 
variable models. 


III. ROBUSTNESS OF GENERALIZED BELL 
INEQUALITIES 

In this Section we characterize as a function of the system 
size the robustness of the above generalized Bell inequalities 
against perturbations of the projections identified in Eqs. & 
This is indeed a fundamental question to address when inves¬ 
tigating the possibility to violate a Bell-like inequality using a 
macroscopic quantum system. In fact, any experimental test 
of this kind will have to cope with non-idealities in the iden¬ 
tification of the most suitable projections to implement. Our 
study is thus aimed at determining the tolerance allowed by 
the Bell-like tests addressed above against imperfect measure¬ 
ment settings. 

In order to attack this problem, we consider the equiva¬ 
lence between an improperly arranged measurement setting 
[i.e., projections of the ideal generalized Bell state performed 
on states that differ from those in Eqs. (|6j, but still being mu¬ 
tually orthogonal] and the performance of ideal measurements 



FIG. 2. (Color online) Histogram of the probabilities P(I) to obtain 
the values I for the maximally entangled state given in Eq. with 
/ = 1 for a perturbation e ~ 0.233. The red line 1 = 3 indicates the 
border between classical (I < 3) and non-classical correlations (/ > 
3). 

on a state that is locally rotated by an unknown amount and 
around an unspecified direction. We thus perform a stability 
analysis in the regime of small perturbations by independently 
rotating each subsystem by an angle e around a random axis. 
The rotations are performed by applying the bi-local unitary 
transformation 

U = U A ®U B = e' E ^ A ® e'^B , ( 7 ) 

where (p = A,B) is a ri-dimensional Hermitian operator 
whose entries are randomly chosen from a uniform dis¬ 

tribution of values within the range — 1 < mu < i. 

We have performed an extensive numerical analysis of the 
value of the Bell expressions I and I,] within the framework 
discussed above. In order to obtain the probability distribu¬ 
tions P(B) for B = I,Id we have accumulated statistics by 
means of N = 10 5 random rotations for each value of e and 
d considered in our study. A typical example of the results 
of this analysis is reported, for B = I, in Fig. [2] the distribu¬ 
tion P(I) is a Gaussian-like distribution with a negative skew. 
In order to provide a quantitative characterization of the ro¬ 
bustness of the Bell inequalities, we focus on the violation 
probability 

P B = f P(B’)dB (8) 
Jb c 

where B c = 2,3 for B. B 1 = Id-,1, respectively. Such quantity 
is the area of the distribution P(B) at the rightmost side of the 
LHV value B c . 

We have studied the behavior of the the probability P B 
against the dimension of the problem at hand: in Fig. [3] (top 
panel) we plot P B as a function of the spin quantum number / 
(d = 2/ + 1), for various choices of e and for both I and Id- We 
first note that for / = 1/2 the results for the two Bell expres¬ 
sions coincide. On the contrary, from / = 1 onwards, Pj d > Pj, 
namely violation of the Bell inequality for Id, which includes 
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FIG. 3. (Color online) Upper panel: Violation probabil¬ 
ity ‘I’b as a function of the spin quantum number l for £ = 
0.12,0.18,0.23,0.29,0.41 (green triangles, blue squares, black di¬ 
amonds, magenta circles, red filled circles). Solid and dashed curves 
show the fit 0 for 7 and I d , respectively. Lower panel: estimate from 
Eq. d 10k of the threshold lj as a function of £, for different threshold 
probabilities (see legend). 


number of terms growing with the dimension of the system, is 
more resistant to local perturbations than violation for 7. It is 
clear from the top panel of Fig.[3]that for all values of £ that we 
have considered and for a given minimum violation threshold 
one can find a threshold value I b (e, Tf,) such that <n 
for all 1 > lg. In order to quantify /^, we fit the data of Fig. [3] 
(top panel) with the function 


1 

2 


1 — erf 


Me) A ’ 


(9) 


characterized by the mean value /#(e) and the standard devia¬ 
tion A/;(e). This allows us to estimate V K as 

IKzXb) « / B (e)+A B (e)err 1 (l -2S§), (10) 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 


£ 

FIG. 4. (Color online) Maximal value / max (e) for different spin 
quantum number /. Red horizontal line indicates the violation bound 
/ = 3. 


with l B (e) « 0.17e-°- 61 and A B (e) « 0.137e^ 138 for 5=7; 
l B (e) « 0.057e~ 2 - 52 and A B (e) « 0.484e^ a50 for B = I d . In 
Fig. [3] (bottom panel) we report Ij for several values of 'Bf. 
These results clearly indicate that even for small error e, 
the violation probability becomes rapidly negligible when in¬ 
creasing l. 

Besides studying the probability to violate a given inequal¬ 
ity, it is interesting to consider whether or not one is actually 
able to violate a Bell inequality in a ^-dimensional system, 
for a given e and a finite sampling. This is equivalent to esti¬ 
mate the quality of the statistical sampling of the tails of the 
Gaussian-like distribution P(B). In order to do so, we adopt 
a pragmatic approach according to which we consider it pos¬ 
sible to violate inequality for the Bell expression B if, with a 
sampling of N = 10 s events, a single instance of violation is 
found. We then extract the maximum occurred value 77 max of 
the Bell expression and study it against 8 for increasing sys¬ 
tem dimensions. Fig.[4]shows the results for the case of B = 7, 
highlighting the existence of a value e) above which the corre¬ 
sponding inequality is no longer violated. Although the quan¬ 
titative analysis is not reported here, we have found that an 
analogous critical value 8^ exists for the case of 77 = I d , too. 

The values of both 8/ and at given / are reported in 
Fig. [5] which also displays the comparison with the fitting 
function a/l b , thus showing a power-law decay of the error 
threshold with the size of the system, with b sa —1.13 for 7 
and b « —0.96 for l d . Consistently with our findings on the 
probability to violate a given inequality, we find that I d is more 
robust than 7 as the corresponding critical value E* f decays 
with l more slowly than . 
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We complete this Section by studying the stability of the 
generalized Bell inequalities when global perturbations are 
used instead of the bi-local ones considered so far. By repeat¬ 
ing the same analysis illustrated above under the assumption 
that the random perturbation now acts on both the spins of 
our system, thus possibly deteriorating their shared correla¬ 
tions, we find the results reported by the green lines in Fig. [5] 
A global rotation of the state of the system results in a more 
pronounced decay of the threshold value of e that, however, 
does not affect significantly the scaling with the size of the 
system. From the same power law fit a / l b as above we obtain 
b « — 1.39 for I and b ~ — 1.57 for Ij. For random global ro¬ 
tations, it can be shown analytically (see appendix^ that the 
asymptotic (in d ) decay of of the violation probability fP/ rf is 
at least exponentially fast. 


IV. RANDOM MEASUREMENT OPERATORS 

In this section, we consider the violation of generalized Bell 
inequalities by a bipartite spin-r/ system prepared in an un¬ 
known state. In order to make our study more specific and 
concrete, we first consider a spin system prepared in a maxi¬ 
mally entangled state but defined over random z-axes. In such 
a case, it is no longer possible to identify the form of the opti¬ 
mal measurement operators As and /is, and the measurement 
settings given in Eqs. ([6) would be unsuitable for the success¬ 
ful falsification of a given generalized Bell inequality. In order 
to address this point quantitatively, we resort to the observa¬ 
tion that the picture given above is fully equivalent to that of 
a spin system prepared in the ideal generalized Bell state and 
subjected to unknown local measurement settings. 

We thus perform measurements along random directions 
and check if this procedure allows to detect a violation with a 
statistically relevant probability. This procedure is not equiv¬ 
alent to the analysis we have done in the previous Section, as 
the direction of the projections entailed by the measurement 
steps needed for the construction of the Bell function are now 
fully random and not mutually orthogonal. Every measure¬ 
ment operator is rotated individually. 

In Fig. [6] we report the resulting distributions of the val¬ 
ues taken by I for various dimensions of the system. The 
distributions can be well fitted by Gaussian functions spec¬ 
ified by mean value 7 and standard deviation o decreasing 
with l. Specifically, we find the best fit values 7 l~ 0 79 and 
o oc r 1 - 13 . This implies that the probability of measuring a 
violation is as small as fP « 4 • 10~ 7 for / = 1, and is below 
machine precision for any l > 1. This is an interesting result, 
as it entails that, for systems with four or more local states, it 
is virtually impossible to violate a generalized Bell inequal¬ 
ity, even when operating over a maximally entangled state but 
with random direction. 

Given that random search fails to detect violation in gen¬ 
eral, we adopt a more sophisticated strategy and optimize the 
observables to be measured, employing the Nelder-Mead Sim¬ 
plex algorithm ED. which is a general-purpose gradient-free 
algorithm of vast use in complex multi-variable searches. We 
repeat the optimization from different randomly chosen initial 



FIG. 5. Log-log plot of critical e^(/) with corresponding fits (form 
a/l b ) for B = I (solid red with squares) and B = 7/ (dashed red with 
circles). The solid green line with squares denotes the global critical 
Ej g for global rotations for inequality /, the dashed green line with 
circles is ejjj for inequality !<]. 



FIG. 6 . (Color online) Histograms P(I) of the Bell inequality I for 
Bell state |<f> + ) and random measurement operators A\, A 2 , B\, and 
62 . f° r various spin numbers /. 


sets of coefficients for the matrix representation of the observ¬ 
ables As and Bs and report in Fig. [7] the maximum achieved 
value of I and 7/. As it can be seen, in both cases a violation 
of the Bell inequalities is found for systems of up to l = 4.5 
and there is indication (as it should be) that for all system 
sizes it should be possible to find a set of observables capa¬ 
ble to detect the violation. However, the computational time 
of one function evaluation scales as 0(f) for inequality I and 
as 0(f) for 7/ fl9l . In combination with an increasing num¬ 
ber of function evaluations needed for increasing /, this once 
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FIG. 7. (Color online) Optimized values / opt (upper panel) and /° pt 
(lower panel) for different spin quantum numbers l (blue curve) and 
maximally entangled state |d> + ). In both panels the bottom (green) 
line corresponds to the border between classical and non-classical 
correlations. 

more practically prevents the effectiveness of such protocol 
for increasing system sizes l20l . 

As a final step in our study, we now check how the violation 
of the generalized Bell inequalities behaves for random two- 
spin states. We considered two cases: 



FIG. 8. (Color online) Optimal values for Bell expressions / (up¬ 
per panel) and Id (lower panel), for random entangled states (blue 
curve) and random separable states (crosses). Green line: bound for 
classical correlations, Triangles: values for / and Ij for the maxi¬ 
mally entangled state and optimal measurement setting as reported 
in Fig. [T] Note: Blue points are always the maximal value of the five 
random states as described above. 


demonstrating once more the increasing difficulties in finding 
a violation for increasing system size. 


1. Random entangled states whose vectors have elements 
randomly chosen from a uniform distribution of values 
within the range [—10,10] and then properly normal¬ 
ized. 

2. Product states of two random d-dimensional states 
(built up as described above). 

For each value of l we consider five different random entan¬ 
gled states and two random product states and store the max¬ 
imally violating instance for both cases. Typical results for 
maximally violating random entangled and a random product 
state are shown in Fig. [8] where we observe that violations are 
possible for both I and Id and arbitrary values of / (blue cir¬ 
cles). This of course is only true for random entangled states, 
as for random product states (red crosses) a violation is impos¬ 
sible. Nevertheless, the critical bounds can always be reached 
by classical random states. Note that once again, finding the 
proper setting to violate the Bell inequalities get harder for 
bigger system sizes: the cases for l > 3.5 in Fig. [8] display 
a lower violation with respect to the other ones. However, 
this is due only to the fact that we stop the optimization af¬ 
ter one week of computational time, which in these cases is 
clearly not enough to achieve the maximal possible violation. 


V. CONCLUSIONS 

We have analyzed quantitatively the violation of Bell in¬ 
equalities formulated for general d-dimensional systems. Our 
approach, which is oriented towards the provision of infor¬ 
mation relevant to experimentalists, was based on the careful 
evaluation of the impact that imperfections in both the state 
resource and the measurement settings could have on actu¬ 
ally violating one of such generalized non-locality tests. We 
have found that both bi-local and global imperfections in the 
arrangements for the measurement settings needed to test a 
given inequality affect the quality of the corresponding Bell 
test in a system size-dependent fashion. Non-locality in bipar¬ 
tite systems spanning large Hilbert spaces is difficult to infer, 
even when managing maximally entangled states of two spin- 
/ particles. Interestingly, we have identified size-dependent 
thresholds on the amount of noise that the system is able to 
tolerate before rendering the chance of violating a generalized 
Bell inequality negligibly small. We have also highlighted the 
influences due to the statistical sampling used to check the fal¬ 
sification of the Bell tests addressed in our work: while a uni¬ 
form random sampling of the measurement sampling makes it 
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unlikely to single out an instance of violation even for qutrit 
systems, a more sophisticated statistical sampling reveals that 
it should be possible to find such instances regardless of the 
size of the local subsystems, although with increasing efforts 
being required. 

We believe that this work helps addressing the inherent 
complexity of violating Bell inequalities above and beyond 
the usual studies used to assess their robustness, and pro¬ 
vides useful information for any experimental endeavor aimed 
at testing non-local features in the state of systems spanning 
large Hilbert spaces, a situation that is quickly becoming more 
realistic and addressable. Furthermore, our work provides for 
a transparent interpretation of the emergence of classicality in 
terms of local realism at the macroscopic level: without re¬ 
sorting to decoherence mechanisms, the size-dependent noise 
thresholds we identified above describe the practical impossi¬ 
bility to detect non-classical behavior in large enough systems 
via Bell inequalities violations. In a word, in the unavoidable 
presence of noise, large systems will behave for all practical 
purposes classically under a Bell test. 

Numerous avenues are left open by the questions addressed 
in this work, the most natural being whether the use of non- 
maximally entangled states leads to a less demanding infer¬ 
ence of non-locality in large-sized systems ED. In this con¬ 
text, it would be interesting to explore the robustness of the 
so-called embezzling states m. i.e. of a family of partially 
entangled states which exhibit maximal violation of Bell in¬ 
equalities, under the random errors addressed herein. 
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Appendix A: Asymptotic decay of tPj d 

Id, as defined in Eq. Q, is a function of the (pure) state t|/ of 
the bipartite system, hence a function on the 2 dr — 1 dimen¬ 
sional real unit sphere . It will be shown that this func¬ 

tion is subject to the “concentration of measure” phenomenon 
H23I . meaning that, if \j/ is chosen at random from the uni¬ 
form distribution on S 2( p_i> then the probability that |/ £ /(\|/)| 
is larger than an arbitrarily small positive number is exponen¬ 
tially small with d\ notably, 

g{feWI>E)<2exp (AD 


To show this it is convenient to rewrite Id in the form: 

4/(¥) = RdM - S d {\|/) , 

= E 5>(*) ll^ ,J Vn 2 , 

l<i,j<2k=0 

S d (w) = E X>(k) 11^1 vll 2 , (A 2 ) 

1<ij<2A=0 

where s = [d/ 2] , c(k ) = 1 —2 k/(d— 1). For i,j = 1,2 
and — d +1 < k < d — 1, is the projector onto the d- 
dimensional subspace wherein A,- — Bj = n(i,j,k), where : 

n(l,l,k) = k, n(l,2,k) = —k, 

n(2,l,k) = —k— \ , n(2,2,k) = k. 

The average of I d (y) with respect to the uniform, normalized 
measure on S 2d 2 _i is 0, because projectors which share the 
same superscript i , j in R d and in S d are unitarily equivalent 
and the uniform measure is unitarily invariant. Then Levy’s 
lemma states that l23l 

P{\I d (v)\>£} < ’ (A3) 

where id is the Lipschitz constant of I d , to be presently esti¬ 
mated. First we write: 

Mw) - j/(40i < \RdM-Rdm + \s d (y)-s d m, 

(A4) 

and we note that, for an arbitrary projector P and arbitrary 
normalized states t|/, <f> 

IHiVl! 2 - ||^I| 2 I < 2||p(i|/—0)||. 

It follows that: 

Mv)-4/(401 

<2 £ £c(*)(||^ , ' ,/) (v-4>)ll + ||^£t(v-40ID- 

l<i.;<2A=0 

( A5) 

In the above sum, projectors that share the superscript 
(i, j) with different k are mutually orthogonal, so P^'’ 1 + 
P^'Jl | is again a projector and 

I l^i‘ J) (¥ - 0) 11 + \\P-kl t (¥ - 4>) II < v / 2||R[ ! ' ;) (¥-0)ll; 

using this in ( |A5[ > we get 

l4/(¥)-4/(401 < 2 V 2 £ E c Wll4' J) (v-0)ll 

l<i,j<2k=0 

y/ 2 /- V /2 

<2Vi e E^w ) EiiA j, (v-«ii 2 

l<i,;<2 \A-=0 / \k=0 / 

< 8^2(1 +s/3) 1/2 ||¥-0l! ■ (A6) 






8 


In the last line, we used that R ‘ l: with the same i,j and 
different k are orthogonal. This inequality shows that i d < 
8\/2(l +S/3) 1 / 2 . Using this in Levy’s lemma yields the an¬ 
nounced estimate. 

A sharper bound is obtained by using the following version 
of Levy’s lemma: 


1/2. While we have no theoretical estimate for M d , numerical 
results suggest that, at least with the observables we have cho¬ 
sen, the distribution of I d is even. This would imply M d = 0, 
and then: 

r P{\I d \>e} < exp(-3c/ 2 E 2 /(128(6Ac/))) . 


nVd~M d | > e} < exp(-r/ 2 e 2 /(2£ 2 )) , (A7) 

where the median M d is defined such that Prob { I,] < M d } = 


While even this stricter bound becomes relevant only for 
dimensions d > 10, bound (All proves rigorously that the vi¬ 


olation probability fP/ rf , upper bounded by the right-hand side 
of (|A1|> for e = 2, decays at least exponentially fast with d. 
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